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UNIVERSAL LOWER BOUNDS FOR POTENTIAL ENERGY OF 

SPHERICAL CODES 

P. G. BOYVALENKOV + , P. D. DRAGNEV D. P. HARDIN*, E. B. SAFE*, 

AND M. M. STOYANOVA** 


Abstract. We derive and investigate lower bounds for the potential energy of finite 
spherical point sets (spherical codes). Our bounds are optimal in the following sense 
- they cannot be improved by employing polynomials of the same or lower degrees in 
the Delsarte-Yudin method. However, improvements are sometimes possible and we 
provide a necessary and sufficient condition for the existence of such better bounds. All 
our bounds can be obtained in a unified manner that does not depend on the potential 
function, provided the potential is given by an absolutely monotone function of the inner 
product between pairs of points, and this is the reason for us to call them universal. 
We also establish a criterion for a given code of dimension n and cardinality N not to 
be LP-universally optimal, e.g. we show that two codes conjectured by Ballinger et al 
to be universally optimal are not LP-universally optimal. 


1. Introduction 

Minimal energy configurations, maximal codes, and spherical designs have wide rang¬ 
ing applications in various helds of science, such as crystallography, nanotechnology, ma¬ 
terial science, information theory, wireless communications, etc. In this article we shall 
derive lower bounds on the potential energy of such configurations via a unified method 
working for a large class of potential interaction functions. A fundamental connection 
between our lower bounds and the classical Delsarte-Goethals-Seidel bounds on spherical 
designs and Levenshtein’s bounds on maximal codes is presented. For a fixed dimension 
and code cardinality the Delsarte-Goethals-Seidel bounds serve to localize the analysis 
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and then, as illustrated in Figure 2, the zeros of the Levenshtein optimal polynomials for 
maximal codes determine the optimal polynomials for a large class of potentials. 

Following Levenshtein’s terminology (see m) we call the lower bounds that we obtain 
universal. This choice of terms is also consistent with its use by Cohn and Kumar in 
their study [TT] of universally optimal energy configurations, since our bounds likewise 
work for all absolutely monotone potential functions of the inner product. Furthermore, 
our lower bounds are attained for all sharp configurations as defined in |14j . 

Let denote the unit sphere in W^. We refer to a finite set C C as a spherical 
code and, for a given (extended real-valued) function h{t) : [—1,1] —>• [0, +oo], we define 
the h-energy of a spherical code C by 

EiC-,h):= h{{x,y)), 

x,y£C,x^y 

where (rr, y) denotes the inner product of x and y. Note that for x,y £ we have 
\x-y\‘^ = 2- 2{x,y). 

A commonly arising problem is to minimize the potential energy provided the cardi¬ 
nality I Cl of C is fixed; that is, to determine 

£{n, iV; h) := inf{K(C; h) : |C| = iV, C C 5^“^} 

the minimum possible h-energy of a spherical code of cardinality N (see |2n[ 131)1. Al¬ 
though the theorems in Section 2 hold for general potentials h we will be especially con¬ 
cerned with functions h{t) that are absolutely monotone {absolutely strictly monotone), 
that is h^^\t) > 0, z = 0,1,... {h^^\t) > 0, i = 0,1,...). Some examples of absolutely 
monotone potentials include the Riesz a-potential h{t) = [2(1 — a > 0, and in 

particular the Newton potential (when a = n — 2); the Gauss potential h{t) = the 

Korevaar potential h{t) = (1 -|- 0 < r < 1. Although the logarithmic 

potential h{t) = —(1/2) ln(l — t) is not positive on [—1,0], all its derivatives are positive 
and the results in this article apply to this potential as well. The situation is similar for 
the Fejes-Toth potential h{t) = —[2(1 — 0 < a < 2, which includes the impor¬ 

tant particular case in discrete geometry of a = 1, namely of hnding configurations that 
maximize the sum of all mutual distances. 

A general technique (referred to here as the Delsarte- Yudin method) for obtaining lower 
bounds for the h-energy of arbitrary spherical codes was developed by Yudin 133 using 
Delsarte’s linear programming method [miisiEi] and was further applied by Kolushov 
and Yudin [52], Andreev |T|, and Cohn and Kumar M- These bounds depend on the 
choice of polynomials satisfying certain constraints. Here we provide explicit solutions 
to Delsarte’s linear program based upon Levenshtein’s work on maximal codes [26| and 
I2ZI, which allows us to establish universal lower bounds on potential energy for a large 
class of potential functions h. 

In Section 2 we describe in a unified manner results from Delsarte, Goethals and Seidel 
|18j and Levenshtein |25l|26l|27| that are instrumental in defining our bounds. Theorems 
2.3 and 2.6 explain the importance of special type quadrature rules in determining lower 
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bounds on energy and investigation of their optimality. Theorem |3.1| is one of the main 
results in this paper. It gives lower bounds which are optimal in the following sense - 
they cannot be improved by polynomials of the same or lower degree that satisfy the 
standard linear programming constraints specified in Theorem |2.2[ On the other hand, 


the bounds of Theorem 3.1 can be further improved in some cases and Theorem 4.1 gives 
necessary and sufficient conditions for existence of such improvements via the so-called 
test functions, which were first introduced and investigated for analysis of the Levenshtein 
bounds for maximal codes in 1996 by Boyvalenkov, Danev and Bumova m- We derive 
a quantitative version of m Theorem 5.2] in Theorem |4.10[ which provides a criterion 
for disproving that certain codes are LP-universally optimal. As an application we prove 
that the two codes conjectured to be universally optimal in [2], are not LP-universally 
optimal, namely their universal optimality may not be established by an ad-hoc approach 
similar to the 600-cell approach given in [laiis]. 


2. Linear programming framework and 1/A1-quadrature rules 


2.1. Gegenbauer polynomials and the Delsarte-Yudin linear programming 
framework. For hxed dimension n, the Gegenbauer polynomials [ 33 ] are defined by 

= 1, = t and the three-term recurrence relation 

{i + n — = {2i + n — 2)tpj^^\t) — for i > 1. 

We note that {Pi^\t)} are orthogonal in [—1,1] with respect to the weight (1 — t2j(«--3)/2 

and that = 1. In standard Jacobi polynomial notation (see |33[ Chapter 4]), we 

have that 


( 1 ) 


, , p{in-3)/2,{n-3)/2), . 

* p((n-3)/2,(n-3)/2)^^^- 


Denote the space of real polynomials of degree at most k by Any f € can 
be uniquely expanded in terms of the Gegenbauer polynomials as f{t) = Yli=o 
The coefficients fi given by 


h 


fll f(t)-Pi(t)(^ - t2)(’^-3)/2 dt 
f-i \pi'"\t)Y {I -dt' 




play an important role in linear programming theorems. 

Let {Yki{x) : f = 1, 2,..., be an orthonormal basis of the space Harm(A:) of homo¬ 
geneous harmonic polynomials in n variables of degree k restricted to where 

^ ,,, /n-|-/c — 3\2A:-|-n — 2 (n + k — 1\ [n + k — Z\ 

r, := d.m HarmW = (^ „-2 j - k -=( n-1 n-1 j 

and orthonormality is with respect to integration over the sphere utilizing an, the nor¬ 
malized (n — l)-dimensional Hausdorff measure restricted to The functions 
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i are known as spherical harmonics of degree k. The Gegenbauer poly¬ 

nomials and spherical harmonics are related through the well-known Addition Formula 
(see [23]): 


( 2 ) 


— '^yke{x)Yke{y) = P]^\{x,y)), x,y e 

Th ^ ' 

^ D _1 


that is, the Gegenbauer polynomial is, up to a normalization, the kernel for the 

orthogonal projection onto Harm(A;). 

If / is a function integrable on [—1,1] with respect to the weight function (1 —12^(«-3)/2 
and y is any fixed point on then the following relation (a partial case of the Funk- 

Hecke formula, see |291 Theorem 6]) holds: 


/gn-l 


f{{x,y))dan{x) =jn / 


where 


In ■■ = 


-1 


r(f) 


V5r(afi)' 

If C = {xi,... ,XAr} is a spherical code of N points on then it follows from Q 

that: 

N rfe AT rfe / AT \ 2 

(3) yk£ixi)Ykeixj) = — X] '^Ykiixi) > 0. 


ru ^ ^ ^ ^ rk 

i,j=l '^e=li,j=l '' £=1 \i=l 


We define the A:-th moment of C by 


N 

Mk{C) := Y.Pt^\{x,,Xj)). 


hi=i 


From (§, we have Mk{C) = 0 if and only if Y{xi) = 0 for all spherical harmonics 
Y G Harm(/c). If Mk{C) = 0 for 1 < fe < r, then C is called a spherical t- design. 
Equivalently, C is a spherical r-design if and only if 



p{x)dan{x) 


1 


x&C 


{an is the normalized (n — l)-dimensional Hausdorff measure) holds for all polynomials 
p{x) = p(xi, X 2 ,..., Xn) of degree at most r. The set 

(4) X(C) :={fcGN: Mfc(G) =0}, 

is called the index set of C. Hence, C is a spherical r-design if and only if {1, 2,..., r} C 
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Suppose / : [—1,1] —)■ M is of the form 

OO 

(5) f{t) = ^ fkPl"'\t), fk>0 for all k>l, 


k=0 


where we remark that /(I) = YlT=o fk < oo- Since \Pl"'\t)\ < 1, it follows that the 
right-hand side of (j^ converges uniformly on [—1,1]. We then obtain the following 
relations which form the basis for many packing and energy bounds for spherical codes 
C = of cardinality N (see PiniEUES]): 


( 6 ) 


N 

i^(C;/)= j;/((x„x,))-/(l)Ar 

OO N 

k=0 

OO 

= Y.fkMk{C)-f{l)N 

k=0 

> foN^ - f{l)N. 

Since Mk{C) = 0 for k = 1,... ,r when C is a r-design, the following result immedi¬ 
ately follows from (|^. 

Theorem 2.1 (Delsarte, Goethals, Seidel [H]). Suppose C is a spherical r-design on 
and f{t) is a polynomial of degree at most r such that f{t) > 0 on [—1,1] and 
fo = In Xli /(t)(l - dt > 0. Then 


( 7 ) 


|C|> 


m 

fo 


Maximizing the right hand side of Q over polynomials satisfying the above hypotheses, 
Delsarte, Goethals, and Seidel [18] obtain a lower bound on 

B{n,T) := min{|C| : C C is a spherical r-design} 

Specifically, they show 

^n + k — 2 


(8) B{n,T) > D{n,T) ;= < 


n — 1 

n -|- fc — 1 
n — 1 




n-\-k — 2 
n — 1 


if r = 2A: — 1, 


if r = 2k. 


We refer to D{n,T) as the Dels arte-Goethals-Seidel bound for spherical r-designs. 
Another application of ([^ is Yudin’s lower bound on energy. 
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Theorem 2.2 (Yudin [35]). Suppose f : [—1,1] —)■ M zs of the form (§ with fk>0 for 
all k > 1. Then, for N >2 

£{n,N;f)>foN^-f{l)N. 

Consequently, if h : [—1,1] —)■ [0,oo] satisfies h{t) > f{t), t G [—1,1], we have 
(9) £{n,N-h)>foN^-f{l)N. 

Furthermore, C is an optimal (energy minimizing) code for h and equality holds in 
if and only if both of the following conditions hold: 

(a) f{t) = h{t) for all t £ {{x,y) : x ^ y, x,y £ (7}; 

(b) for all k > 1, either fk = 0or Mk{C) = 0. 

For a given h : [—1,1] —)• [0, oo], we denote by An^h the set of functions f < h satisfying 
the conditions Q. Recall that for such /, the coefficient sequence (/o, /i,. . .) G £i. The 
problem of maximizing the lower bound /o-Y^ — /(1)Y arising in Theorem |2.2| can then 
be expressed in terms of an infinite linear program; 


( 10 ) 


maximize F{fo, /i, ...):= Y ^/o(Y - 1) - ^ , 

OO 

subject to fkPk^\t) < h{t),t £ [—1,1] and fk > 0, for all A: > 1. 
k=0 


In the following we shall consider the above linear program restricted to a subspace A 
(usually finite-dimensional) of the linear space C'([—1,1]) of real-valued functions contin¬ 
uous on [—1,1], For such a A, we define 


(11) W(n,Y,A;h):= sup N^{fo - f{l)/N). 

/eAnA„,h 


In general, it can be a difficult problem to find the value of W(n, Y, A; h). We consider 
sufficient conditions that allow us to solve for W{n, N, A; h). In particular, we explicitly 
find the solutions of the truncated linear program (10) and thus find 0 when A = Vk, 
for all k < T{n, Y), for some T{n, N) (as defined in equation (18) below). In the particular 
case when m = T{n,N) we derive the universal lower bound (ULB) for potential energy 
of spherical codes. 


2.2. 1/Y-Quadrature rules and lower bounds for energy. We refer to a finite 
sequence of ordered pairs {(oz, /Oi)}f=i as a 1 /N-quadrature rule if —1 < oi < 02 < • • • < 
Ofc < 1, and Pi > 0 for i = 1, 2,..., fe, and say that {(oi, is exact for a subspace 

AcC([-l,l])if 




/(I) 


k 

i=l 


( 12 ) 


Y 
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for all / G A. 

Theorem 2.3. Let {{ai, be a l/N-quadrature rule that is exact for a subspace 

AcC([-l,l]). 


(a) If f £ An An^h, then 


S{n,N]h) > N‘^'^pif{ai). 


1=1 


(b) We have 


(13) 


W{n,N,A-h) <N^Y.Pih{ai). 


2=1 


If there is some / G A n An^h such that f{ai) = h{ai) for i = 1,..., k, then equality holds 
in (13) which yields the universal lower bound 


(14) 


£{n,N-,h) > N‘^'^pih{ai). 


2=1 


Proof. If / G A, then (12) holds and so, from Theorem 2.2, we obtain 


£{n, N; h) > N\fo - f{l)/N) = 

2=1 

showing that (a) holds. 

For (b), using (12), we obtain 

W(n,iV,A;h)= sup N\fo-f{l)/N) 

k k 

= sup 

Clearly equality holds if there is some / G A n An^h such that /(oi) = h{ai) for i = 
1,... ,k. □ 

As we next describe, a spherical code C = {xi,..., xn} C provides a quadrature 
rule that is exact on the subspace 


oo 


Ac:=\m = fo+ Y m Y i/^i< 

\ iex(c) iex(c) 

with I{C) as defined in (Q. Let 

{(xj, Xj ): Xi / Xj G C} =: {-1 < ai < a2 < ■ ■ ■ < at < 1}, 
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and let {qi} denote the inner product distribution; i.e., 

{xi,Xj) = ai}\ 

qi . ^2 ) ‘ L, . . . , K. 

If / G Ac, then /; = 0 for all I 0 I{C) (unless I = 0) and equality holds in Q. Hence, 
for such /, we obtain 


(15) 


/o = ^ {E{C-, f) + Nfil)) = M + qifim), 

1=1 


N 


that is, {{ai, is a l/A'-quadrature rule exact for Ac- 


Example 2.4. As an example we consider the 600-cell C consisting of 120 points in 
Each X € C has 12 nearest neighbors forming an icosahedron (the Voronoi cells 
are dodecahedra) and there are 8 inner products — 1 = ai < a2 < ••• <03 < 1 
between distinct points in C. If f{t) < h{t) on [—1,1] and f{ak) = /i(afc) and for 
all Q!fc > —1, then we must also have f'{ak) = h'{ak), resulting in 2 • 7 -|- 1 = 15 
interpolation conditions. If C were a 14-design, then this would suggest we search for 
/ G A 4 /jnA with A = 7^14. However, C is only an 11-design (i.e., Mi 2 {C) / 0), although 
A7 i3(C) = ••• = Mig{C) = 0, so C is almost a 19-design. This suggests we choose 
A to be a 15-dimensional subspace of "Pig H In fact, Cohn and Kumar |14| 

Section 7] show that for any absolutely monotone potential h on [—1,1], there is a unique 
/ £ An,h n A for A := {/ G P 17 : /n = /12 = /13 = 0} that proves the optimality of C. 


Example 2.5. Another example is provided by the so-called sharp configurations m, 
namely configurations with k distinct inner products that are spherical designs of strength 
2A; — 1. In this case A = V 2 k-i and the existence of the 1/A^-quadrature is provided by the 
configuration quadrature (15) and the design property. We shall return to this example 
in the Remark 3.3 following Theorem |3.1[ 


The two examples above cover all currently known universally optimal conhgurations. 
The next theorem provides sufficient conditions for optimality of (14) even in a larger 
subspace. 


Theorem 2.6. Let {(oi ,be a 1/N-quadrature rule that is exact for a subspace 
A C C'([—1,1]) and such that equality holds in ( [T^ . Suppose A' = A0 span : j G 
/} for some index set I C N. If := ^ -|- Yli=i ^ d for j G I, then 

k 

W(n, iV, A'; h) = W{n, N, A; h) = N^Y^ pih{ai). 

i=l 


Proof. Suppose f{t) G H A'. Then we may write the decomposition of / as 

fit) = 9{t) + 
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for some g € A and fj > 0, for j G I. Note that /o = go, since 0 0/. Furthermore, since 
the quadrature rule {{(Xi, pi)}^^^ is exact for <7 G A, we have 

/o-/(l)iV-i = go-f{l)N-^ = ^-^ + Y,p^g{a,)-{g{l) + Y,fl]N-^ 


2 = 1 




I /(“<) -I - I E/d 


2=1 

k 








jW/ 


i=l 

k 


j&I 


N 

(n) 


2 = 1 
k 




2=1 




2 = 1 


where, for the last inequality, we used /(t) G /j and > 0. 


□ 


2.3. Levenshtein bounds for spherical codes. Let 

A(n, s) := max{|C|: C C {x,y) < s, x ^ y ^ C} 

denote the maximal possible cardinality of a spherical code on of prescribed maximal 
inner product s. 

For a,b £ {0,1} and i > 1, let denote the greatest zero of the adjacent Jacobi 

polynomial Pl ^ ^ (t) and also define t^’ = —1. For r G N, let Ir denote the 

interval 


I-r := < 


1,1 1,0 

‘'k 


1,0 1,1 
''fc ’ ‘'k 


, if T = 2/ — 1, 


if T = 2k, 


The collection of intervals is well defined from the interlacing properties 

see 1211 Lemmas 5.29, 5.30]. Note also that it partitions Z = [—1,1) into countably many 

subintervals with non-overlapping interiors. 

For every s £ Z^, using linear programming bounds for special polynomials 
of degree r (see [271 Equations (5.81) and (5.82)]), Levenshtein proved that (see p7l 
Equation (6.12)]) 


(16) A{n, s) < < 


r _ ^k+n—3\ [2k+n—3 

L 2 k-i{n,s)-{ -. 


if s G X 2 fc-i 


r (y. — (k+n-2\ r2fc+n-l (l+^K-Pfe -Pfe+}(s))l - r T 1 

L 2 k{n,s)-{ , ifsGX 2 U 
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For every fixed dimension n each bound Lrin, s) is smooth with respect to s. The function 

L2k-iin,s), ifsGX2fc_i, 

L{n, s) = 


[ L 2 kin,s), if s Gl 2 k 

is continuous in s. The connection between the Delsarte-Goethals-Seidel bound Q and 
the Levenshtein bounds (flG]) is given by the equalities 


(17) 


L2k-2{n,tl’\) 

L2k-i{n,tl’^) 


and the ends of the intervals Ir 


= L 2 k-iin,tl’\) = D{n,2k - 1), 
= L 2 k{n,t]f) = D{n,2k). 



Figure 1. The Levenshtein function L(4, s) on 1 < A; < 6. 


2.4. Levenshtein’s 1/A^-quadrature rule. Levenshtein’s method for obtaining his 
bounds on cardinality of maximal spherical codes utilizes orthogonal polynomials theory 
and Gauss-type quadrature rules that we now briefly review. The location of the cardi¬ 
nality N relative to the Delsarte-Goethals-Seidel numbers D{n,T) is an important step 
in determining our universal lower bounds. From the properties of the bounds D(n, r) 
and LT-{n, s) (see (§, (0) we derive that for every fixed dimension n and cardinality N 
there is unique 

(18) T:=T{n,N) such that N ^ {D{n,T), D{n,T + 1)]. 

For the so found r define k := \and let ak = s he the unique solution of 

(19) N = Lr{n,s), s G Xt-. 
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Then as described by Levenshtein in I2ZI Section 5] (see also [26l E] ) there exist uniquely 
determined quadrature nodes and nonnegative weights 

(20) - 1 < ai < • • • < Ofc < 1, pi, ... G M+, i = l,...,k 

such that the Radau/Lobato 1/Ai-quadrature (see [16], |5]) holds 


( 21 ) 


/o — 


m 

N 


+ ^ Pi/(«i), for all / G Vr- 


2=1 


When T = 2/c — 2 is even, then ai = —1 and (21) is Lobato quadrature. The numbers 
dj, i = 2,... ,k, are the roots of the equation 


( 22 ) 




Pk-i{t)Pk-2{oik) - Pk-i{ak)Pk-2{t) = 0 , 


where Pi{t) = pi ^ ^ '[t). When r = 2/c — 1 is odd, then ai > —1 and (21) becomes 

Radau quadrature. The numbers ai, i = 1..., k, are the roots of the equation 


(23) Pk{t)Pk-i{o:k) - Pkio:k)Pk-iit) = 0, 

where Pi{t) = P- ^ ^ \t). In fact, {oi} are roots of the Levenshtein’s polynomials 

/r"’’“'')(t) (see [271 Equations (5.81) and (5.82)]). 

The dynamical behavior of the quadrature nodes {a^} is the following. When N G 
{D{n, 2k —2), D{n, 2k—1)) then ai = —1 and the quadrature (21) is Lobato. The solution 


ttfc of (19) belongs to the interval and all {ai }^^2 strictly increase with N. 

We have that 

1 = |ai| > \a2\ > \ak\ > laal > |ciA:-i| > • • • . 


ii.i 


At the transition point N = D{n,2k — 1), ai = —1 and ak = The equation (22) 

, which implies that 

1 = |ai| > \ot2\ = Wk\ > |«3| = lufc-il > • • • . 


becomes T)[”J]^(t) = 0, which implies that 


As N increases from D{n, 2k —1) to D{n, 2k), ak strictly increases from to as do 
the rest of the nodes {ai}^~l. In particular, ai > —1 and (21) defines Radau quadrature 
and 


R.o 


1 > |q:i| > \ak\ > \ct2\ — |<afc_i| > • • • . 

More details on the nodes {ai} can be found in [121 Appendix], [THl Corollary 3.9], and 
[3 Section 2.6]. 


3. Universal lower bounds 


3.1. Optimal polynomials for lower bounds. The optimal polynomials of degrees 
one and two to be applied in Theorem 2.2 can be found by direct computations and 
manipulations with the corresponding derivatives. These polynomials suggest a general 
form of polynomials which are optimal in the following sense - they give lower bounds 
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which cannot be improved by utilizing other polynomials of the same or lower degree in 
Theorem 12.21 


Our choice of polynomials for Theorem 2.2 can be viewed as extension of the ideas of 
Levenshtein [261 [27] who uses suitable quadrature formulas (Subsection 2.4) to explain 
the bounds (16) and their optimality in the same sense as above. This similarity should 
not seem unusual - the maximal code problem is infinite version of the Riesz energy 
problem. In fact, Cohn and Kumar [T4| use similar idea to deal with the universally 


optimal configurations. Thus, our paper can be viewed as natural extension of the works 
[261 EH [HI . Recall that given a fixed dimension n and a code cardinality N we can 
associate r = T{n,N) and s G Xt such that Lr{n, s) = N (see (18) and ([I^). Depending 
on the parity of r we distinguish two cases: 

Case (i): t = 2k —2 and Uk = s G • Then f{t) := Hermite 

interpolation polynomial of degree 2k — 2 defined by (recall that ai = — 1 in this case) 

(24) /(-I) = /i(-l), f{ai) = h{ai), /'(a*) = h'{ai), i = 2,...,k. 


1,1 1,0 
''fc-i’ '-fc 


Case (ii): t = 2k — 1 and ak = s G 
interpolation polynomial of degree 2k — 1 defined by 

(25) f{ai) = h{ai), /'(a*) = h'{ai), i = l,2,...,k; 


Then f{t) := is the Hermite 


In the notation of Cohn-Kumar’s paper m p. 110], our polynomials are 

(26) /*„,„,(*) = (*-»)/'’“!(()). 


3.2. Main theorem. The equations ( [24| ) and ( [2^ define a Hermite’s interpolation prob¬ 
lem for f{t) to intersect and touch the graph of the potential function h{t) (see |22[ The¬ 
orems 2 and 3], [HI Section 5]). This implies as in [HI Sections 3 and 5] that / G An^h 
and we could use f{t) for bounding S{n,N; h) from below. Observe that the nodes ( [20[ ) 
are independent of the potential function /i, hence we call our bound on £{n,N]h) a 
universal lower bound (ULB). 

Next, we state our main theorem. We note that here is the first time when we impose 
the condition that the potential function h{t) is absolutely monotone and that none of 
the preceding results have required this property. 


Theorem 3.1. Let n, N he fixed and h{t) he an absolute monotone potential. Suppose 
that T = T{n,N) is as in (18), and choose k = Associate the quadrature nodes 

and weights ai and pi, i = 1,... ,k, as in (21). Then 


(27) 


k 

£{n, iV; h) > Rr{n, N; h) := ^ pih{ai) 

i=l 


Moreover, the polynomials defined by (i), respectively by (ii), provide the unique optimal 
solution of the linear program 0 for the subspace K = Vr and consequently. 


(28) 


W{n,N,Vr-,h) = Rr{n,N;h). 
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Remark 3.2. The optimality of the Hermite interpolants (26) is analogous to the op¬ 


timality of the Levenshtein polynomials (proved first by Sidelnikov [32j), and 

emphasizes the universality of our bound. 


Remark 3.3. As noted in Example 2.5 
quadrature. Moreover, the k inner products coincide with {a,}. 


the sharp configurations (see [TTj) define 1/N- 

Consequently, the 


bounds (27) are attained by all sharp configurations. 


Proof of Theorem 3.1. We first consider the odd case (ii), that is r = 2^ — 1. The 
conditions in (ii) define Hermite interpolation at the points Oj, i = 1, 2,..., fc, and give 
a unique polynomial / of degree 2A: — 1 with positive leading coefficient. The absolute 
monotonicity of h{t) implies that f{t) < h{t). 

Next we derive that / satisfies the condition ([^ as well. From ( |23[ ) we have that 
the quadrature nodes {ai,... ,ak} are zeros of the polynomial Pk{t) + cPk-i{t), where 

{Pi} are the Jacobi orthogonal polynomials {P^ ' ^ }. From the interlacing proper¬ 

ties of the orthogonal polynomials we obtain that the constant c = —Pk{s)/Pk-i{s) is 
non-negative. Indeed, the largest roots of the Jacobi polynomials of Pk-i satisfy 
tl’-i < ifc-i (see [26]). Since the last but largest root of Pk is smaller than (by 
the interlacing property), we obtain that the ratio Pk{t)/Pk-i{t) doesn’t change sign in 
f^’°). Moreover, from [71 Lemma 3.1.3 (a)] (see also [H Lemma 1.5.8]) 


Pkitl’^i) n + 2k-3 


> 0 , 


n + 2k-l 

hence c > 0. Utilizing the approach of m Sections 3 and 5] we conclude that the Hermite 
interpolant / has non-negative Gegenbauer expansion. Therefore, / G /j. 


We now use (21) to derive the universal bound of /. We have 

k k k 

+ ^(/oAf-/(!)) = 


/o — 


/(I) 


N 


i=l 


2=1 


2=1 


which means that S{n, N; h) > J2i=i Pih{oci) = i? 2 fc-i(^) h). 


Furthermore, for any polynomial u = Y^i=o ^ '’^iPi^^{t) £ ^n,/i of degree at most 2A; — 1 
we have 

k k 

(29) NifoN - /(I)) = N^Y1 ^ 




2 = 1 


2 = 1 


i.e. N{uoN — u{l)) < R2k-i{n, N;h) and u{t) does not improve (27). 

Should equality hold in (29) for some u G n 7^2fc-i) we observe that u{ai) = h{ai) 
for f = 1,2,...,A:. Additionally, the condition u{t) < h{t) implies that u'{ai) = h'{ai) 
for all Qfj G (—1,1). Hence, u satisfies the Hermite interpolation data (25), and by the 
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uniqueness of the Hermite interpolant, u = f. Therefore, / is the unique optimal solution 
to the linear programming problem (10) in the class An^h fl'P 2 fc-i and (28) holds. 

In the even case (i) we proceed analogously, where we only modify the proof of the 
non-negativity of the Gegenbauer expansion. In this case we utilize m Lemma 10]. □ 



Figure 2. The optimal polynomials (Hermite interpolants), that provide 
the ULB for Gauss, Korevaar, and Newton potentials (in ascending order), 
along with the corresponding Levenshtein polynomial for n = 4, = 24 


3.3. Discussion and examples. The bounds (27) are easy for computation and in¬ 
vestigation. Moreover, the approach by which they were derived doesn’t depend on the 
potential function and in this sense they are universal. This universality is illustrated 
in Figure 2, where we consider n = 4, = 24 and plot the Gauss, Korevaar, and New¬ 

ton potential functions, together with the corresponding optimal Hermite interpolants of 
degree r = 5, that solve the linear program (10) in the class V^. We also overlay the 
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Levenshtein polynomial whose zeros are the solutions of (23), where s satisfies 

L5(4, s) = 24. These zeros of the Levenshtein polynomial also serve as quadrature nodes 
for the universal lower bound (27) and as Hermite interpolation nodes for the optimal 
LP polynomials. 


In [2] the authors have done an extensive experimental investigation of energy-minimizing 
point configurations, in particular they provide the computational minimizers for the 
Newton potential energy = [2(1 — when n = 1,2,..., 32 and N = 

1,2,, 64. Table 1 compares the Newton energy from [2] and our universal lower bound 
(ULB) when n = 4 and N = 5, 6,..., 64. 


N 

Harmonic 

Energy 

ULB Bound 

% 


N 

Harmonic 

Energy 

ULB Bound 

% 


N 

Harmonic 

Energy 

ULB Bound 

% 

5 

4.00 

4.00 

0.00 


25 

182.99 

182.38 

0.34 


45 

664.48 

663.00 

0.22 

6 

6.50 

6.42 

1.28 


26 

199.69 

199.00 

0.35 


46 

697.26 

695.40 

0.27 

7 

9.50 

9.42 

0.88 


27 

217.15 

216.38 

0.36 


47 

730.75 

728.60 

0.29 

8 

13.00 

13.00 

0.00 


28 

235.40 

234.50 

0.38 


48 

764.59 

762.60 

0.26 

9 

17.50 

17.33 

0.95 


29 

254.38 

253.38 

0.39 


49 

799.70 

797.40 

0.29 

10 

22.50 

22.33 

0.74 


30 

274.19 

273.00 

0.43 


50 

835.12 

833.00 

0.25 

11 

28.21 

28.00 

0.74 


31 

294.79 

293.51 

0.43 


51 

871.98 

869.40 

0.30 

12 

34.42 

34.33 

0.26 


32 

315.99 

314.80 

0.38 


52 

909.19 

906.60 

0.28 

13 

41.60 

41.33 

0.64 


33 

337.79 

336.86 

0.28 


53 

947.15 

944.60 

0.27 

14 

49.26 

49.00 

0.53 


34 

360.52 

359.70 

0.23 


54 

985.88 

983.40 

0.25 

15 

57.62 

57.48 

0.24 


35 

384.54 

383.31 

0.32 


55 

1025.76 

1023.00 

0.27 

16 

66.95 

66.67 

0.42 


36 

409.07 

407.70 

0.33 


56 

1066.62 

1063.53 

0.29 

17 

76.98 

76.56 

0.54 


37 

434.19 

432.86 

0.31 


57 

1108.17 

1104.88 

0.30 

18 

87.62 

87.17 

0.51 


38 

460.28 

458.80 

0.32 


58 

1150.43 

1147.05 

0.29 

19 

98.95 

98.48 

0.48 


39 

487.25 

485.51 

0.36 


59 

1193.38 

1190.03 

0.28 

20 

110.80 

110.50 

0.27 


40 

514.90 

513.00 

0.37 


60 

1236.91 

1233.83 

0.25 

21 

123.74 

123.37 

0.30 


41 

543.16 

541.40 

0.32 


61 

1281.38 

1278.45 

0.23 

22 

137.52 

137.00 

0.38 


42 

572.16 

570.60 

0.27 


62 

1326.59 

1323.88 

0.20 

23 

152.04 

151.38 

0.44 


43 

601.93 

600.60 

0.22 


63 

1373.09 

1370.13 

0.22 

24 

167.00 

166.50 

0.30 


44 

632.73 

631.40 

0.21 


64 

1420.59 

1417.20 

0.24 


Table 1. Newtonian (harmonic) energy comparison (see [2]) with ULB 
for n = 4, N" = 5,..., 64. 


Utilizing the same Newton energy-minimizing conhgurations provided in [2] in Table 
2 we compare our universal lower bound (ULB) with the Gauss potential {h{t) = e^*“^) 
energies of these configurations, which in general provide upper bounds on the minimal 
Gauss energy for the same choice of n = 4 and N = 5,6,..., 64. We note that the error 
dramatically improves, which is to be expected, as the Hermite interpolants of analytic 
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potential functions are excellent approximants. Observe that for N = 5 and N = 8 the 
bounds are exact. Both cases are universally optimal. 


N 

Gaussian 

Energy 

ULB Bound 

% 


N 

Gaussian 

Energy 

ULB Bound 

% 


N 

Gaussian 

Energy 

ULB 

Bound 

% 

5 

0.82085 

0.82085 

0.0000 


25 

54.83402 

54.81419 

0.0362 


45 

195.4712 

195.46 

0.0042 

6 

1.51674 

1.469024 

3.1460 


26 

59.8395 

59.7986 

0.0684 


46 

204.7676 

204.76 

0.0049 

7 

2.351357 

2.303011 

2.0561 


27 

65.02733 

64.99832 

0.0446 


47 

214.2834 

214.27 

0.0075 

8 

3.3213094 

3.321309 

0.0000 


28 

70.43742 

70.41329 

0.0343 


48 

223.994 

223.99 

0.0007 

9 

4.6742772 

4.614371 

1.2816 


29 

76.06871 

76.0435 

0.0332 


49 

233.9421 

233.93 

0.0040 

10 

6.1625802 

6.123668 

0.6314 


30 

81.9183 

81.88889 

0.0359 


50 

244.0939 

244.09 

0.0022 

11 

7.9137359 

7.85 

0.8517 


31 

87.99142 

87.95307 

0.0436 


51 

254.4665 

254.46 

0.0028 

12 

9.8040902 

9.780806 

0.2375 


32 

94.26767 

94.2326 

0.0372 


52 

265.0585 

265.05 

0.0049 

13 

11.975434 

11.92615 

0.4116 


33 

100.75 

100.7275 

0.0223 


53 

275.8551 

275.85 

0.0030 

14 

14.353614 

14.28178 

0.5005 


34 

107.4465 

107.4377 

0.0082 


54 

286.8694 

286.86 

0.0020 

15 

16.90261 

16.88487 

0.1049 


35 

114.3862 

114.3632 

0.0202 


55 

298.1012 

298.1 

0.0019 

16 

19.742184 

19.70346 

0.1962 


36 

121.5266 

121.504 

0.0186 


56 

309.5522 

309.54 

0.0030 

17 

22.795437 

22.73703 

0.2562 


37 

128.874 

128.86 

0.0109 


57 

321.2188 

321.21 

0.0041 

18 

26.046099 

25.98526 

0.2336 


38 

136.4529 

136.4314 

0.0158 


58 

333.0979 

333.08 

0.0043 

19 

29.510614 

29.44794 

0.2124 


39 

144.244 

144.218 

0.0180 


59 

345.1882 

345.18 

0.0033 

20 

33.161221 

33.12489 

0.1096 


40 

152.2451 

152.2199 

0.0165 


60 

357.497 

357.49 

0.0033 

21 

37.051623 

37.03121 

0.0551 


41 

160.4628 

160.4379 

0.0155 


61 

370.0202 

370.01 

0.0030 

22 

137.52 

137.00 

0.3753 


42 

168.8894 

168.8713 

0.0107 


62 

382.7551 

382.75 

0.0019 

23 

41.177514 

41.15351 

0.0583 


43 

177.5346 

177.5199 

0.0083 


63 

395.7039 

395.7 

0.0004 

24 

45.537431 

45.49154 

0.1008 


44 

186.3928 

186.3839 

0.0048 


64 

408.8804 

408.87 

0.0021 


Table 2. Gauss energy of the harmonic optimal configurations (as pro¬ 
vided in [2]) compared with ULB for n = 4, = 5,..., 64. 


As a consequence of the proof of Theorem 3.1 we describe the explicit LP solutions 
for m < T{n, N) in the next corollary. 


Corollary 3.4. The linear program (LP) can be solved for any m < T{n,N) and the 
solution in the class Vm H An^h is given by the Hermite interpolants at the Levenshtein 
nodes determined by N = Lm{n, s). 


Example 3.5. Here we present the suboptimal LP solutions for n = 4 and N = 24. 
In this case T{n,N) = 5. For m = 1,... ,5 we find the intersection of iV = 24 with 
Li(4, s),..., L5(4, s). The corresponding suboptimal solutions as Gegenbauer expansions 
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(up to three digits) are: 

hit) = .499Po(i) + .229Pi(t) 

/ 2 (t) = .581Po(t) + .305^1 (t)+ 0.093P2(t) 

hit) = .658Po(t) + .395Pi(t) + .183P2(t) + 0.069P3(t) 

hit) = .69Po(t) + .43Pi(t) + .23P2(t) + .10P3(t) + 0.027P4(t) 

hit) = .71Po(t) + .46Pi(t) + .26P2(t) + .13P3(t) + 0.05P4(t) + 0.01P5(t). 



Figure 3. Suboptimal LP solutions for n = 4 and N = 24. 

A natural question is whether linear programming bounds can be improved if we 
consider polynomials of higher than r(n, iV) degree. The next section investigates this 
topic. As one would expect from our results thus far presented, the analogy with the 
situation for maximal spherical codes is quite close. 


4. Necessary and sufficient conditions for optimality of the universal 

LOWER bounds 


4.1. Test functions. Let n and N be fixed, r = T(n, A^) and Lrin^s) = N be as in 


(18) and (19), and j be a positive integer. We introduce the following functions in n and 


s = Ofc: 


1 ^ 

Qj(n, s) ■■= — + '^ piPj'^\ai) for s £ h 
i=l 


( 30 ) 
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It follows that Qj{n, s) = 0 for 1 < j < t and s € Ir (since this is the coefficient /o = 0 
in the Gegenbauer expansion of Thus the functions Qj{n, s) are not interesting 

for these cases and so we assume below that j > t + 1 when s € Ir- 

The next theorem shows that the functions Qj{n, s) give necessary and sufficient con¬ 
ditions for existence of improving polynomials of higher degrees. 


Theorem 4.1. The bounds (21) can be improved by a polynomial from An^h of degree at 


least T + 1 if and only if Qj(nfs) < 0 for some j > r + 1. Furthermore, if h is strictly 
absolutely monotone and Qj{n, s) < 0 for some j >t + 1, then (21) can be improved by 
a polynomial from An^h of degree exactly j. 

Proof. We give a proof for t = 2k — 1. 


(Necessity) The necessity follows from Theorem 2.6 for I = {2k, 2A: -|- 1,... }. 

(Sufficiency) Conversely, assume that h is strictly absolutely monotone and suppose 
that Qj{n, s) < 0 for some j > 2k. 

We shall improve the bound p7|) by using the polynomial 


f{t) = ePj'^\t) + g{t), 

where e > 0 and g{t) G V2k-i will be properly chosen. Denote h{t) := h{t) — eP^^\t) 
and select e such that h{t)^^\t) > 0 on [—1,1] for all i = 0,1,... ,j. Observe, that for 
this choice of e the function h{t) is absolutely monotone. The polynomial g{t) is chosen 
as the Hermite interpolant of h at the nodes {ai}, i.e. 

g{ai) = h{ai), g'{ai) = h!{ai), i = 1,2,... ,k. 

Since h{t) is an absolutely monotone function, we infer as in Theorem 3.1 that g G A^ 
implying that / G An,h- 

Let g{t) = Yl‘e=~o^ Note that /o = go and /(I) = ^(1) + e. We next prove 

that the bound given by f{t) is better that R 2 k-iin, N; h). To this end, we multiply by 
Pi and sum up the first interpolation equalities: 

k k k 


'^Pihioii) = '^Pih{ai) - piP'f‘\ai 


2 = 1 


2=1 


Since 


'^Pig{ai) = go 


2=1 


2=1 


5 ( 1 ) 

N 


by (21) and 




(n) 


2=1 


(ai) = Qj{n,s) - — 
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by the definition of the test functions (30), we obtain 

9(1) 


90 


N 


R 2 k-i{n,N-,h) e ^ / x 
-- + --eQ,{n,s) 


which is equivalent to 

N{Ngo - {g{l) + e)) = R 2 k-i{n, N- h) - eN^Q,{n, s). 

Therefore N{NfQ - /(!)) = R 2 k-i{n,N;h) - eN'^Qj{n,s) > R 2 k-i{n, N;h), i.e. 
polynomial f{t) gives better bound indeed. We also obtained a new bound 

(31) W{n, N; h) > R 2 k-i{n, N; h) - eN^Qj{n, s). 


the 


□ 


Theorem 4.1 provides a sufficient condition for solving the infinite linear program (10) 


Corollary 4.2. If Qj{n,s) > 0 for all j > T{n,N), then {f) solves the linear 

program (10). 


4.2. Investigation of the test functions. The test functions (30) coincide with the 


functions with the same name which were introduced and investigated in 1996 by Boy- 
valenkov, Danev and Bumova m- More details and all proofs are given in the disserta¬ 
tions [7] and [H]. We cite some results from [HEiis] with only reformulations for energy 
bounds. 


Theorem 4.3 ([7j, |8], [H] ). The bounds Rr{n,N] h) cannot he improved by using poly¬ 
nomials of degrees r + 1 and t 2. 


Set A:i(n) := y/n — 2 and let k 2 {n) > 9 be such that 

4n < k 2 {n)‘^ — Ak 2 {n) -|- 5 -|- y/k 2 {nY — 8 k 2 {n)^ — 6 ^ 2 ( 71 )^ -|- 24 fc 2 (n) -|- 25. 

Then we have the following theorems. 

Theorem 4.4. a) [TJ Theorem 3.5.15], [El Theorem 3.4.12] If n > 3 and k > ki{n), then 
all bounds R 2 k{n, N; h) corresponding to s in the open interval l 2 k can be improved by 
polynomials of degree 2 k -\-3. 

b) [71 Theorem 3.5.9], [H Theorem 3.4.14] If n > 3 and k > k 2 {n), then all bounds 
R 2 k-i{n, N] h) corresponding to s in the open interval l 2 k-i can be improved by polyno¬ 
mials of degree 2k + 3. 


Theorem 4.5. a) If n > 3 and k > ki(n), then 

(32) £{n, N] h) > R 2 k-i{n, N; h) - eN‘^Q 2 k+?,{n, s). 

for every N G (D(n,2k — 1), D{n,2k)) where e is chosen as in Theorem 4-1 
b) If n > 3 and k > k 2 {n), then 

(33) S{n, N; h) > i? 2 fc(n, iV; h) - eN‘^Q 2 k+^{n, s). 


for every N G {D{n, 2k), D{n, 2k + 1)) where e is chosen as in Theorem f-l 
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Proof. This follows from (31) and the fact that Theorem 4.4 is based on the inequality 
Q 2 k+ 3 in, s) < 0 which holds true for the mentioned values of n and r. □ 


Another application of Theorem |4.4| concerns the sharp configurations. Recall that a 
sharp configuration is a maximal spherical (n, L 2 k-iin, s), s) =(dimension, cardinality, 
maximal cosine) code; i.e. a code that attains the odd Levenshtein bound L 2 k-i{n, s) 
(cf. |26]). In fact, the next corollary is implicit in 1121 and follows from the main result 
of [28] as well. 

Corollary 4.6. For any fixed dimension n > 3 only finitely many sharp configurations 
are possible. 


Proof. Theorem |4.4| implies that in every fixed dimension n > 3 every Levenshtein 
bound L 2 k-i{n,s) can be improved in the whole open interval provided k is 

large enough. The remaining end points correspond to tight spherical designs, which 
means (among many other things) that /c < 6 [3l 0]. This leaves only finitely many 
possible intervals F 2 k-i where the Levenshtein bound L 2 k-i{n, s) can be attained. Every 
such interval contains finitely many s, corresponding to cardinalities N, which completes 
the proof. □ 

We complete the subsection with the following conjecture, based on the above results 
and numerous investigations of the test functions as related to maximal spherical codes. 

Conjecture 4.7. If Qj{n, s) > 0 for j = T{n, N) + 3 and T{n, N) + 4, then Qj{n, s) > 0 
for all j > T(n, N). 


4.3. Test functions and LP universality. We now apply the test functions to the 
study of universal configurations. 

Definition 4.8. A spherical code C C of cardinality \C\ = N Is called LP- 

universally optimal if 

E{C; h) = >V(n, iV, "P; /i), for all absolutely monotone /i. 


where V is the subspace of polynomials. 


Remark 4.9. Observe that from (§ and 0 one infers that LP-universally optimal 
codes are in fact universally optimal. If the conjecture in Ballinger et al [2] is true, then 
Theorem 4.11 implies that the converse does not hold. 


We derive a criterion for positivity of test functions of large enough j that can be 
used for proving that certain spherical codes of given dimension n and cardinality N are 
not LP-universally optimal. We utilize {n,N) to denot^ codes C C M"" with cardinality 
\C\ = N. As examples, the cases {n,N) = (10,40), (14,64) and (15,128) are analyzed. 


^We note that [2] uses (N, n) notation instead. 
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Sharp estimations for Gegenbauer polynomials can be derived from m (see also 
In |19l Theorem 1] the following inequality is given 

(34) max + 

te[-i,i] TT 

where {pj{t)} are the orthonormal Jacobi polynomials with weight w{t) = (1 —t)“(l+t)^. 


Utilizing a = f3 = to get Gegenbauer polynomials and the normalization = 1, 

we rewrite (34) as 

(35) 


(n) 


Wl< 


r(^) 


2»^-2e(4 + (n-3)\/2)j! 


(1 — t2)(n-2)/4 y 7 r( 2 j + n — 2) (j + n — 3)! ’ 


where r(x) is the Gamma function |34] . Note that for every fixed n > 3 and t G (—1,1) 
the right-hand side of (35) is strictly monotone decreasing in j. 

Let k, ai, a 2 and pi be as in Theorem 3.1 Denote by jo{n,N) the smallest degree 
j > T{n, N) such that the right hand side of (35) is less than when 


(36) 


t = 


if > —1, 

if ai = —1 and pi < ^, 


or less than when t = 02 if «i = —1 and pi = 1/N. 


Theorem 4.10. Let n > 3, N > 2, and let k, ai, 02 o,nd pi be as in Theorem 3.1. Then 
Qj{n, ak) > 0 for all j > jo{n, N). 

Proof. As the comments on the dynamical behavior of the quadrature nodes {a^} at the 
end of Section 2 indicate, we have |ai| > \ai\ for i = 2,..., /c and in the case oi = —1, 
we further have |a 2 | > \oii\ for i = 3,..., fc. 

We first consider the case |ai| < 1; i.e., ai > — 1. If j > jo{n, N) then we have 


(37) 


1 


k 

E 

Z=1 






1 


N 


1 


N 


1 


iV - 1 


= 0 


(we used p* = — 1 following from (12) for f{t) = 1). The case ai = —1 and 

Pi < l/A" is handled similarly using (35) as suggested by the second line of (36). 

For the final special case oi = — 1 and pi = 1/N it is clear (cf. [I2]) that Qj{n, s) = 0 


for odd j. The case of even j follows similarly as above using the facts that Pj"'\—1) = 1 
and that |aj| < |a 2 | for i = 3,..., /c. □ 


Theorem 4.10 gives a useful tool for disproving LP-universal optimality. For given n 
and N and numerics suggesting that Gorollary |4.2| may hold one finds explicit jo{n,N) 
and calculates the remaining test functions Qj{n, s) for every j G {r(n, N) + 3, T{n, N) + 
4,... ,jo{n, N) — 1}. This will be applied in the next subsection for some codes from [2]. 
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4.4. Examples. Table 3 lists the first twenty test functions for some interesting config¬ 
urations. We utilize {n,N) to denote codes C C MT" with cardinality \C\ = N. 


j 

(4, 24) 

(10,40) 

(14, 64) 

(15,128) 

(7,182) 

(4,120) 

0 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

4 

0 

0.021943574 

0.013744273 

0.000659722 

0 

0 

5 

0 

0.043584477 

0.023867606 

0.012122396 

0 

0 

6 

0.085714286 

0.024962302 

0.015879248 

0.010927837 

0 

0 

7 

0.16 

0.015883951 

0.012369147 

0.005957261 

0 

0 

8 

-0.024 

0.026086948 

0.015845575 

0.006751842 

0.022598277 

0 

9 

-0.02048 

0.02824122 

0.016679926 

0.008493915 

0.011864096 

0 

10 

0.064232727 

0.024663991 

0.015516168 

0.00811866 

-0.00835109 

0 

11 

0.036864 

0.024338487 

0.015376208 

0.007630277 

0.003071311 

0 

12 

0.059833108 

0.024442076 

0.01558101 

0.007746238 

0.009459538 

0.053050398 

13 

0.06340608 

0.024976926 

0.015644873 

0.007809405 

0.0065461 

0.066587396 

14 

0.054456422 

0.025919671 

0.015734138 

0.007817465 

0.005369545 

-0.046646712 

15 

-0.003869491 

0.02498472 

0.015637274 

0.007865499 

0.006137772 

-0.018428319 

16 

0.008598724 

0.024214119 

0.015521057 

0.007815602 

0.005268455 

0.020868837 

17 

0.091970863 

0.025123445 

0.01562458 

0.007761374 

0.005134928 

-0.000422871 

18 

0.049262707 

0.025449746 

0.015694798 

0.007812225 

0.004722806 

0.012656294 

19 

0.035330484 

0.024905002 

0.015617497 

0.00784714 

0.003857119 

0.006371173 

20 

0.048230925 

0.024837415 

0.015589583 

0.00781076 

0.007863772 

0.011244953 


Table 3. Test functions for some special {n,N) spherical codes. 


Judging by the behavior of the test functions the linear programming method will 
provide improvements on our ULB for (4, 24) and (7,182) but it is unlikely to give a 
solution similar to the case with the 600-cell (4,120), where a polynomial in Vn served 
as an exact lower bound. Indeed, that the test functions Qii,Qi 5 and Qn are negative 
provides additional insight on the unique property of the 600-cell as the only universally 
optimal code known that is not a sharp configuration. 

The first configuration (4, 24) is the root system, or the so-called kissing number 
configuration in (see |30) ). which was shown by Cohn, Conway, Elkies, and Kumar 
(see [13]) not to be universal. The negative test functions s) and <59(4, s), s = 02 ~ 
0.4749504897, suggest searching for a polynomial f{t) = Y^^=o /e = /? = 0 

and four touching points of the graphs of f{t) and the potential h{t). We have developed 
a numerical algorithm for handling such situations. For example, if h{t) = 2{i-t) 
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Newton potential, our numerical calculations led to the polynomial 

f{t) = 0.4987 + 0.4852t + 0A535f + 0.5546t^ + 0.9401t'^ + 0.8425t® 

- 0.3305t® - 0.7479t^ + 0.1889t® + 0.37394t^ 

= 0.0073P9^'‘^(t) + 0.0066P8^^^(t) + 0.0659P5^^^(t) + 0.2384Pi'^^(t) + 0.5116Pj'^^(t) 

+ 0.7915^2^'^^ (t) + 0.9236^1^"^^ (t) + 0.7U2P^'^\t). 

The Hermite interpolation points are approximately —0.860297, —0.489872, —0.195724 
and 0.47850. The bound obtained from f{t) is 333.1575, while the universal lower bound 
pT] ) gives ^ 5 ( 4 , 24; 1/(2(1 — t)) = 333 and the energy of the II 4 root system is 334. 
Theoretical and computational aspects of the aforementioned algorithm for improvements 
(when possible) of our ULB and their nature will be discussed elsewhere. 

Theorem 4.11. The spherical codes {n,N) = (10,40), (14,64) and (15,128) are not 
LP-universally optimal. 


Proof. The codes (10,40) and (14,64) were conjectured by Ballinger, Blekherman, Cohn, 
Giansiracusa, Kelly, and Schiirmann in [2] to be universally optimal. It follows from 
Theorem 4.10 and numerical calculations as explained in the end of the last subsection 
that these codes are not LP-universally optimal. Indeed, we have r(10,40) = 3 (so 
ai > —1), jo(10,40) = 10 and the second column in Table 3 shows that this code is 
not LP-universally optimal. Similarly, r(14, 64) = 3, jo(14, 64) = 8, and the inspection 
of the third column of Table 3 suffices. The code (15,128) was not conjectured to be 
universally optimal (but not eliminated) in [2] and we see that it is not LP-universally 
optimal because of r(15,128) = 3, jo(14,64) = 9, and the fourth column in Table 3. □ 
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